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PRESSURE AND HEAT TRANSFER ON B LUNT~CUR VED PLATES WITH 
CONCAVE AND CONVEX SURFACES IN HYPERSONIC FLOW 

By Davis H. Crawford 
Langley Research Center 

SUMMARY 

Pressure and heat transfer have been measured in the presence of both favorable 
and adverse pressure gradients at Mach numbers of 6.8 and 9.6 and Reynolds numbers 
based on length from 0.5 x 10® to 3.0 x 10®. The pressure data were represented by a 
power law relationship, and the empirical constants from this relationship were used to 
calculate the heat transfer by a simplified method. This method allows similar com- 
pressible solutions of the laminar boundary layer with pressure gradient to be directly 
applied without the use of integral or iterative methods. The results of these calcula- 
tions are compared with the experimental values and with the calculated results from a 
more complex laminar similarity method requiring a machine program and iterative pro- 
cedures. Differences between the two theoretical methods are shown to be slight for 
many applications. Applications where the simple method is less accurate are defined. 

INTRODUCTION 

Boundary- layer calculations by similarity methods are applicable with accuracy to 
a wide range of flow conditions where deviations from a calorically perfect gas are small. 
These similar solutions may be applied with good results to a diverse assortment of 
hydrodynamic configurations. Even when the caloric imperfections of the fluids are not 
small, the results of similar solutions are considered conservative when used to estimate 
heat transfer, displacement effects, and momentum thickness. 

Similarity solutions of the boundary layer are ordinarily tabulated as functions of 
quantities in the transformed plane where similarity was achieved. The use of these cal- 
culations involves a thorough understanding of the derivation of the relations as well as 
tedious numerical operations which are best assigned to a high-speed electronic computer. 
Some references have reported these solutions in a form usable in a straightforward 
manner which circumvents the intricate digital procedure. 

One such similarity method by Li and Nagamatsu presented boundary-layer calcu- 
lations for hypersonic flow. (See ref. 1.) This scheme used p °c x n and the hypersonic 



approximations, u e ~ Uoo and 1 - ^ ~ 1, to simplify the boundary -layer equations in 

t 0 

the transformed plane. The assumption of similarity, which allowed the momentum equa- 
tion to be reduced to an ordinary differential equation with the independent variable r), 
resulted in a relation between the Falkner-Skan velocity gradient parameter /3 and the 
pressure gradient parameter n. A greatly simplified method by Bertram and Feller 
(ref. 2) has been advanced for laminar boundary -layer calculations which used this hyper- 
sonic relation between /3 and n with a different Np r = 1 similar compressible solu- 
tion of the boundary layer to correct a suitable flat-plate solution for the effects of pres- 
sure gradient. Bertram shows that this method is inapplicable when n < -1, and implies 
that this method is of doubtful accuracy when the value of n is near -1. 

Another method recently introduced by the present author (ref. 3) provides solutions 
of the boundary -layer equations in closed form. This method which does not involve 
hypersonic approximations is applicable when the pressure distribution is of the form 
p oc ( x + A) n . Assumption of this type of pressure distribution along with a perfect gas 
equation of state and constant entropy at the edge of the boundary layer allowed the 
Np r = 1 compressible solutions of references 4 and 5 to be expressed directly in terms 
of quantities in the real plane. When the hypersonic approximations (t e /t 0 ) — 0 and 
(too/to) — 0 are made, and when A — 0, the relations of reference 3 become identical to 
the relations of reference 2. Comparison of values calculated by these two methods can 
therefore be used to determine the effect of hypersonic approximations for specific flow 
conditions. 

In the present investigation, measured pressure and heat transfer were obtained on 
a two-dimensional model with a hemicylindrically blunted leading edge and one concave 
and one convex curved surface. The main body of the data was taken with end plates, but 
some data were taken without end plates to show the end-plate effect. The model was 
tested at M = 6.8 and 9.6 and at a = -8°, 0°, and 15° for both Mach numbers. The 
free-stream Reynolds numbers based on model length varied from 1 x 10® to 2.8 x 10® 
at M = 6.8 and from 0.6 x 10® to 1 x 10® at M = 9.6. Therefore, laminar flow condi- 
tions occurred for most of the tests. 

This study was initiated to verify the effect of a pressure gradient on the boundary 
layer. This investigation involved situations in which the simple theory of reference 2 
was inadequate since the pressure was poorly represented by the exponential law p oc x n 
with a constant value for n, and the best approximate fit of the data for constant n pro- 
duced values of n near the troublesome value of -1. Also, the local flow velocity was 
often less than hypersonic. The assumption of a local similarity concept, and thus a 
locally determined value of n, allowed the simple theory of reference 2 to be applied only 
on the leading 30 percent of the test area (ref. 3), and in this instance it did not agree 
with the trend of the data although there was agreement with the general level of the data. 
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The slightly more general exponential pressure law used in the method of the 
present author (ref. 3) extended its applicability to better cover the experimental cases 
involved in this investigation. To verify more completely the applicability of this method, 
the calculated heat -transfer results were compared with those computed by a well- 
established more complicated theory which involves a point -by-point machine solution of 
the boundary -layer equations of reference 6 with similarity assumed. 

SYMBOLS 


constant in assumed exponential pressure variation 

velocity of sound 

constant 

drag coefficient of leading edge 

total skin-friction coefficient on a flat plate, -y ^ Cf dx 


C D,£e 

Cf 

Cf 

c w 

C P 

c p 

d 

f 

G 

H 

h 

h>h 


local skin-friction coefficient 


linear viscosity coefficient, 


linear viscosity coefficient, 


P oo^w 

p'too 

Poot' 


specific heat at constant pressure 
diameter of hemicylindrical leading edge 
function of 77 related to stream function 

function of wall temperature and specific heat ratio in laminar -boundary - 
layer growth equation 

,,2 

total enthalpy, h -h -ii— 
local enthalpy, Cpt 

integrals which represent solutions of boundary -layer equations 


i 



Kg local inclination of boundary layer with respect to plate surface in hypersonic 

similarity form, M w d6/dx 

K 1 ,K 2 ,K 3 ,K 4 ,K 5 factors used to correct flat-plate boundary-layer calculations for 

pressure gradient effect 

l total length of flate plate or total length of model 

M Mach number 

m pressure gradient parameter, U e = X m 

Np r Prandtl number 

Ngj- Stanton number 

n constant in assumed exponential pressure variation 

P ratio of local wall pressure to free-stream pressure, p^p^ 

p pressure 


Pb 

P v 

Pw 

R 


S 


t 


U 


u 

X 


pressure evaluated for blunt flat plate at 0° angle of attack 

pressure caused by boundary-layer growth 

pressure evaluated for sharp flat plate at local inclination 


Reynolds number 

enthalpy function, - 1 
He 

temperature 


an 

transformed velocity, u 

a e 


air velocity relative to test surface 


px pa 

transformed longitudinal coordinate, \ X — dx 

JO p o a o 
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Y 

y 

a 

P 

y 

6 


V 

e 


a e r y P 
a o Po 


longitudinal coordinate 
transformed normal coordinate 
normal coordinate 
angle of attack 

pressure gradient parameter, 
ratio of specific heats 
boundary -layer thickness 
boundary -layer displacement thickness 


2m 
m + 1 


function of y, 1 - 


0.0048 

(r - !) 2 


function defined in equation (A3) 
similarity variable, 


, Yi|m + 1 U e X 

Xf 2 v Q 




boundary -layer momentum thickness 

ut n 

linear viscosity coefficient, -jj— 
dynamic viscosity 


kinematic viscosity 


P 

Xcc 

Subscripts: 


density of air 


viscous interaction parameter, 


Moo 3 ^ 


evaluated at edge of boimdary layer 


i 


evaluated at reference position 



evaluated at trailing edge of flat plate 


l 

o evaluated at free-stream stagnation conditions 

tr transformed quantities 

w evaluated at wall conditions 

x evaluated at x distance from leading edge 

°o evaluated at free-stream conditions 

Primes denote differentiation with respect to tj or parameters evaluated by the 
T-prime method. 


THEORY 

Similarity has been achieved on a hypersonic transformed plane in reference 1 by 
assuming p oc x n . The relation between /3 and n shown in this reference has been 
used in reference 2 to evolve a very simple method for determining heat transfer, skin 
friction, and boundary-layer thickness (which has been shown nearly equal to displace- 
ment thickness in hypersonic flow) in hypersonic flow with a pressure gradient. In this 
same reference, however, it is demonstrated that the solutions become imaginary, and 
thus invalid, as n — -1. 

The details of a theoretical method by the present author for calculating character- 
istics of the laminar boundary layer are shown in reference 3. This method involves the 
more general assumption of p <x (x + A) n and does not require the assumption of hyper- 
sonic flow conditions. These assumptions, along with boundary -layer growth in a con- 
stant entropy field and the perfect gas equation of state, are applied to the Np r = 1 
solutions of references 4 and 5 to produce solutions in closed form. Calculations by this 
method are shown in reference 3 which do not become imaginary for a pressure gradient 
described by a value of n less than -1, and for local flow conditions less than hyper- 
sonic. A summary of the theoretical results shown in reference 3 are given in 
appendix A. 

Experimental pressures were used to determine the constants A and n utilized 
in the theory of reference 3. Theoretical estimates of these pressures by a simple 
hypersonic method were tried with varying degrees of success. Explanation of this 
method and the results are shown in appendix B. 
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Subsequent applications of relations given in reference 3 showed that, for small val- 
ues of A and for values of n approaching -1, the calculated parameters were unreal- 
istically sensitive to small changes in the value of A. This application is approaching 
the area for hypersonic flow where the relations of reference 3 reduce to those of refer- 
ence 2, and where those of reference 2 are imaginary. A study of the parameters involved 
shows that when n is negative, the assumption p a x n produces an infinite pressure 

on the stagnation line. Even more unrealistic is the fact that when n approaches - 

3y - 1 

(-7/8 for air), the integral for the longitudinal coordinate transformation used in refer- 
ences 2 to 5 (x = f \ dx) fails to converge. The modified assumption p cc (x + A) n 
\ d 0 Po a o ) 

prevents the infinite pressure predicted at x = 0 for a negative n when A has a posi- 
tive finite value and also prevents the divergence of the integral in the longitudinal coor- 
dinate transformation. The parameter arising from this integral which applies to the 
theory directly is W - J- 1 2L when A = 0 and t e /t 0 = 0, the value of this parameter 
may be expressed as follows: 


m + 1 x dX 
2 X dx 


n + 1 


It is apparent that this relation shows no discontinuities or irregularities in the neighbor- 
hood of n = -7/8. Thus, another consideration must be examined to learn the conse- 
quences of the nonconvergence of this integral when A = 0 and t e /t Q = 0. 

If the pressure is assumed to vary as p cc (x + A) n , then 4^ = n P or 
x dP n dxx+A 

p dx = — J ~K’ ^ a cons ^ an * Pressure and pressure gradient are considered for a sta- 

* + x 

tion xp p -gE- is constant and n is a function of A at this station x^. The effect of 


l* 


a change in the values of A jx^ and n on the shape of the pressure curve when n is 

negative is illustrated in figure 1(a). It is clear that families of curves mutually tangent 

at a point all have the same value of at that point, and that curves of each family 

with increasing curvature have increasing values of A/x^ and increasing negative values 
of n. The effect that p A/x^, and n have upon the heat transfer is given in fig- 
ure 1(b) where the heat -transfer parameter has been calculated for the hypersonic case 
(t e /t 0 = 0), for y = 1.4, and for Npj. = 0.7. The heat transfer is reduced for increasingly 
negative values of p Also, at all values of p the effect of A/x^ upon the heat 
transfer is reduced at large values of A/xp At values of p near zero, the effect of 

A jxi upon the heat transfer is small at all values of A/xp but at values of p 


approaching 


~ 2 r 


(-7/8 for air), the values of calculated heat -transfer parameters 


are very sensitive to small changes in A/x^ when A/x^ is near zero. The rapid 
change in these calculated values is apparently beyond that which should be expected from 
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physical considerations. This effect is the apparent consequence of the previously men- 
tioned singular integral. This characteristic of the theory will be recalled subsequently 
during the comparison of theory and data. 

EXPERIMENTAL APPARATUS 
Tunnel and Test Conditions 

The experimental temperature pressure and visual flow measurements were made 
in the Langley 11-inch hypersonic tunnel in the M = 7 and M = 10 air nozzles. The 
Mach 7 nozzle is a two-dimensional nozzle constructed of Invar since this metal has a 
very low temperature coefficient of expansion. Thus, the temperature strain of the first 
minimum during a test was low and the variation of Mach number with time was less than 
that expected for an identical steel nozzle. The unit Reynolds number for these tests 
was from 5 x 10® per meter to 14 x 10® per meter. The Mach 10 nozzle is a three- 
dimensional nozzle so that the temperature strain near the first minimum has little effect 
on its effective area, and the Mach number does not vary during a run with constant stag- 
nation conditions. The unit Reynolds number per meter for the tests in this nozzle varied 
from 2.5 x 10® to 5 x 10®. The air supplied these nozzles is dried by being pumped 
through a silica gel drier as it is stored at 50 atmospheres pressure before it enters the 
storage tank. During a test the air enters the nozzle after passing through an adjustable 
pressure-regulating valve, an electrical heater with nickel-chromium-alloy heating ele- 
ments, and a settling chamber where the vorticity and contamination of the air are reduced 
by screening. An adjustable second minimum and a diffuser followed by a water-cooled 
heat exchanger downstream of the test region delays the pressure rise in the vacuum tank 
and thus delays the breakdown of the flow in the test nozzle. Many of the details of con- 
struction, calibration, and instrumentation of this intermittent facility are reported in 
references 7 to 10. 

The Reynolds numbers of the tests in the M = 7 nozzle were approximately 
1.0 x 10®, 1.5 x 10®, 2.0 x 10®, and 2.8 x 10® based on free-stream conditions and the 
total length of the curved plate. The Mach number for these tests varied slightly with 
tunnel position, stagnation pressure, and elapsed time from the start of the test but always 
fell in the range 6.6 to 6.9. The Mach number variation with time and with stagnation 
pressure was accounted for in the calculation of the pressure and heat -transfer param- 
eters. The stagnation temperatures fell in the range between 585° K and 655° K but 
varied much less than this during a single test. 

The Reynolds numbers of the tests in the M = 10 nozzle were approximately 
0.5 x 10® and 1.0 X 10® based on free-stream conditions and the total length of the model. 
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The Mach number was between 9.5 and 9.75 and varied with tunnel position and stagnation 
pressure. The stagnation temperature for these tests varied from 820° K to 880° K. 

Models and Instrumentation 

The model profile used in this research program was a cylindrical leading edge 
followed by identical curved surfaces with one convex face and one concave face exposed 
to the airstream. The curvature of these surfaces was determined as shown in figure 2 
from the relation y/B = (x/B)®/®, started with a slope of 7r/3 (x/B = 0.0662), and con- 
tinued until x/B reached a value of 5/8 and the surface inclination was ir/5. This arc 
was the generator of the airfoil's lower surface and was rotated n/18 and translated a 
sufficient distance to join the nose profile to become the airfoil's upper surface. The 
nose profile was a cylinder of 6.35-mm diameter followed by wedge surfaces, 1.6 mm 
long with 7t/18 included angle, tangent with the two airfoil surfaces. The position of this 
wing resting on a flat surface was defined as the zero angle of attack of the wing, the base 
was normal to this flat surface, and the upper surface was extended to meet it. 

In this investigation two models were tested at angles of attack of -8°, 0°, and 15°. 
The positive angle was defined with the concave surface exposed to the flow. The three 
angles which were chosen were such that the concave surface near the trailing edge was 
alined with the stream at a = -8°, the convex surface near the trailing edge was alined 
with the stream at a = 0°, and the convex surface near the leading edge was alined with 
the stream at a = 15°. 

One of the test models was constructed of solid aluminum and contained no instru- 
mentation. This model, shown in figure 3, was sting mounted and was used to study the 
shock formation and surface shear by means of schlieren and oil-flow studies. 

The other model was constructed of 1.27-mm-thick Inconel sheet mounted on a steel 
frame as shown in figure 4. Temperature and pressure measuring stations were provided 
on both the concave and convex surfaces of this model, and removable end plates were 
supplied. The numerous pressure and temperature leads required the model to be side 
mounted. A photograph of the model in place is shown in figure 5. 

The transient skin temperatures were measured by chromel-alumel 30-gage 
(0.254-mm-diam.) thermocouple wire spot-welded to the inner surface of the 1.27-mm- 
thick skin of the model. The two wires approaching each measuring station were shaped 
so that they diverged and then turned toward each other on an axis normal to the air flow. 
The two wires were spot -welded in this position not more than 1.6 mm apart. Thus, the 
output of each thermocouple station represented the average temperature over a lateral 
distance of as much as 1.6 mm but over a longitudinal distance of less than 0.25 mm (the 
diameter of the thermocouple wire). The thermoelectric signals from these stations were 



indicated and recorded on fast-response 18-channel galvanometer oscillographs. The 
cold junctions were located in a junction box in close proximity with a thermometer and 
shielded from rapid changes in ambient temperature by insulation. Calibration of the 
galvanometers was accomplished by impressing measured currents through precision 
resistances in series with the measuring stations and the galvanometers. 

Pressure orifices 1.52 mm in diameter were secured in place with silver solder 
at each pressure -measuring station. Leads no longer than 125 cm were connected with 
the six-cell aneroid-type recording instruments described in reference 7. The test cham- 
ber was evacuated prior to the test to prevent pressure errors from outgassing of the 
lines. The total temperature in the settling chamber was measured by bare wire ther- 
mocouples located at different radial stations and recorded on self-balancing 12 -channel 
recorders. Stagnation pressures were observed on high-accuracy Bourdon tube type of 
pressure gages as well as on one cell of the fast-response instruments previously men- 
tioned. The film record of the total pressure was examined to confirm instant establish- 
ment of the test pressure and to synchronize the timing of the test measurements. 

Temperature Data Reduction 

The heat transfer to this model was determined from the transient temperature 
response of the model skin near the start of the test. This method is explained in detail 
in reference 11. Sample conduction calculations indicate that conduction contributions 
were no more than 1 percent of the measured heating. Therefore, the experimental 
results are presented without conduction corrections. The recovery temperatures used 
to reduce local heating rates to local film coefficients of heat transfer were calculated 
from measured local pressures. Flow conditions at the outside edge were determined 
from an assumed isentropic expansion from the total pressure behind a normal shock. 

The recovery factors used were the square root of the local Prandtl numbers evaluated 
at temperatures calculated from Monaghan's T' method (ref. 12). 

RESULTS AND DISCUSSION 

The nature of the flow about this blunted curved plate may be learned from the 
schlieren pictures of the flow about the uninstrumented sting -mounted model with no end 
plates (fig. 6). The detached shock which wraps itself around the body is nearly sym- 
metric about an axis not alined with the free stream. For the present tests, this shock 
envelope is more nearly alined with the free stream when a = 15°. At this angle of 
attack, the breadth of the shock envelope is a minimum. 

There is little optical evidence that the boundary layer is separated either on the 
convex or on the concave surface of the wing. The schlieren photographs show a weak 
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shock formation initiated as a coalescence of compression waves emanating from the con- 
cave surface of the wing about three diameters behind the leading edge. The expansion 
waves at the base of the vehicle are those expected for an attached boundary layer. 

Two schlieren photographs are presented for flow at a = 15°, M = 6.8, and 
R l = l .3 x 10®. Both photographs are data from the same test at different times. 

Although the shock about the model appears the same in both, tunnel wall separation is 
evident in one and not in the other. Such unsteady flow phenomena have previously been 
reported in situations with separated flow (refs. 11 and 13). When the model with the 
pressure and temperature instrumentation was in place, the tunnel blockage was greater 
than when the visual flow sting-mounted model was in place. Therefore, when the pres- 
sure data from some tests exhibited rapid fluctuations and failed to correlate well, the 
results from these tests were assumed to be faulty because of tunnel interference and 
were discarded. 


End Plates 

Oil -flow photographs in figure 7 show the direction and relative intensity of the sur- 
face shear. The flow on the plate caused by the tip effect is shown clearly in the oil -flow 
photographs. Along the center line where the instrumentation is concentrated, the end 
flow appears to be small, and the use of end plates might be considered unnecessary. 
Nevertheless, in an attempt to obtain better two-dimensional flow over the test surfaces, 
end plates were used while most of the data were obtained. Subsequent oil-flow photo- 
graphs show that weak shocks caused by shock boundary -layer interaction in corner flow 
where the end plates attach to the wing also interfere with two-dimensional flow. 

The heat transfer and pressure data may be examined to show the difference 
accorded by the presence or absence of end plates. Figure 8 was prepared to compare 
these results, which are shown only for M = 6.8. The data at M = 9.6 followed the 
same trend as the lower Mach number data with a little more scatter. On the convex 
surface, the pressure data with end plates fall above the data without end plates. At -8° 
angle of attack, the data are evidently increased about 20 percent near the trailing edge. 

At 0° angle of attack, the magnitude of the effect upon pressure was about the same. At 
15° angle of attack, the pressure data gave evidence of separation on the convex surface. 
This evidence was seen both with and without the end plates. The separation appears 
somewhat more extensive with the end plates. The effect of the end plates on the heat 
transfer to the convex surface is partly obscured by scatter, but appears to be less than 
that on pressure. Most of the end-plate effect on the heat-transfer parameter at a = -8° 
and at a = 0° is in the opposite direction of the effect on the pressure. Thus, the 
parameter Ngj. oo\J^x } co can be expected to correlate the data better for a = -8° and 0°. 
For the 15° angle of attack, however, the parameter used correlated the heat -transfer 
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data very well even though the pressures vary by as much as a factor of 2. Both with and 
without end plates, however, the heat -transfer data for a = 15° indicate a small region 
of laminar separation followed by transition. On the concave surface the effect of the end 
plates on the pressure and the heat-transfer data is not as clear as the effect on the data 
on the convex surface. At a = -8° the end plates apparently cause the data with end 
plates to scatter near the trailing edge. This scatter is possibly caused by transitional 
flow along one row of thermocouples and laminar flow along another row. At a = 0° 
the end plates have more effect on the pressure and less on the heat transfer, whereas 
at a = 15° the effect of the end plates on both the pressure and heat transfer is small. 

There is little question that weak shocks generated by end plates affected part of 
the experimental results by changing the pressure slightly and by affecting the location of 
the start of transition. Their absence, however, led to the signal seen in the oil-flow data 
fostered by the tip vorticity and transmitted obliquely along the surface. Thus, the bene- 
fit of using end plates is unresolved. 


Convex Surface 

The experimental pressures on the convex surface of the curved plate are shown in 
figure 9. End plates were used in this and all subsequent data in this report. The trend 
of these data is well represented by the empirical power law pressure distributions which 
were used to calculate the heat -transfer theory. The experimental and theoretical heat- 
transfer results for the corresponding cases are also shown in figure 9. Most pressure 
data for a given test Mach number and angle of attack appear little affected by the range 
of Reynolds numbers used in the tests, and most of the heat-transfer results are likewise 
little affected by the Reynolds number variation. 

At -8° and 0° angles of attack, no part of the convex surface was at a negative 
inclination to the stream. The base pressure was below stream and there was thus no 
downstream pressure rise to be transmitted forward through the boundary layer. The 
pressure decreases monotonically from the leading edge to the trailing edge. The exper- 
imental heat-transfer results agree in general level and trend with the laminar theory. 

The boundary layer is thus established as laminar and attached. 

At 15° angle of attack, the convex surface was inclined at negative angles 
approaching 15° at stations toward the trailing edge. The possibility of forward trans- 
mittal of a downstream pressure rise was thus presented by this configuration. The 
resulting evidence of separation at M = 6.8 is seen in figure 9(e). The deviation from 
the pressure trend expected for attached flow indicates possible separation over the last 
one-third of the wing's surface. Only a hint of the anticipated reduction in heat transfer 
usually associated with laminar separation is observed, however, and this apparently 
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occurs near the start of separation. The steep increase in heat transfer after separation 
indicates the possibility of transition on the less stable separated boundary layer. 

The experimental pressure results for a = 15° determined at M = 9.6 for the 
convex surface in figure 9(f) do not show separation with the same degree of certainty 
as was shown for M = 6.8. However, the reduced heating on the final one-third of the 
wing surface is an indication of the thickening of the boundary layer just prior to separa- 
tion, if not the actual separation. The pressure never reaches free stream as it did at 
the lower Mach number due to the blast -wave effect of the leading edge. The higher 
pressure on the concave side of the wing caused the higher base pressure which was 
transmitted forward on the convex surface. The sharp increase in heat transfer seen 
just ahead of the trailing edge at the greater Reynolds number again signifies the possible 
start of transition on a separated boundary layer. 


A detailed examination of the theoretical and experimental heat-transfer results 
for figures 9(a) to 9(e) shows good agreement between the theory of reference 3 and the 
elaborate machine calculation. The data fall from only slightly below the theory to 
approximately 30 percent below the theory. The same comparison between the theoret- 
ical values calculated by the two methods and between the theoretical values and the data 
was expected for a = 15° on the convex surface at M = 9.6. The greater scatter in 
the pressure data for this experimental case increased the uncertainty in determining 
the power law pressure function for use in the heat -transfer theory. As seen in fig- 
ure 9(f), a straight line appeared to fit the data within the tolerance suggested by the 
scatter. For this data fit, A = 0, and the negative value of n determined was only 


slightly less than 0.875, the value of 


2y 


for air. The resulting theoretical values 


3y - 1 

of the heat -transfer parameter determined by the two methods from this set of conditions 
failed to show the same agreement and to exhibit the correlation with the data expected 
from previous results. As a first attempt to understand the lack of agreement between 
the two theoretical methods, a different power law variation of pressure was determined 
with an arbitrary value of 1 for A, and with the curve constrained to pass through two 
points near the ends of the previous p versus x function. As may be observed from 
figure 9(f), this pressure curve is very slightly different from that with A = 0 and 
agrees with the pressure data as well as the first trial. The negative value of n speci- 
fied is slightly greater than 0.875 rather than slightly less. The theoretical heat-transfer 
results determined by the theory of the present author for this slight perturbation in 
faired pressures was increased by approximately 50 percent and moved from a position 
below the data to one above the data. The values determined by the method of reference 6 
were only slightly affected by this trivial change in the pressure distribution and were 
thus bracketed by the two theoretical results realized from the simple theory. 
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The extreme sensitivity of the values calculated by the method of the present author 
(ref. 3) to small changes in the value of A when A is small is discussed in the section 
entitled "Theory. M Notice that the limiting values of p/p^ as x approaches 0 for 
these two assumed power law pressure functions are °° and 39, whereas the correct 
stagnation line value is 119. When a third power law pressure function was devised which 
was constrained to pass through the same end points and to predict the correct pressure 
on the stagnation line, it followed a path between those of the first two pressure functions. 
The heat -transfer values calculated with this third power law pressure function are in 
good agreement with the results of the complete machine solution of the boundary -layer 
equations of reference 6, and the theories agree well with the data. This agreement may 
be taken as an indication that, as A tends toward zero and n tends toward -2y/(3y-l), 
the closed-form solution of reference 3 will give reliable results if the assumed power 
law pressure function is constrained to predict the correct pressure on the stagnation 
line. 


Concave Surface 

The pressure distribution and heat transfer on the concave surface of the blunted 
curved plate is seen in figure 10. The effect of the changing Reynolds number has little 
effect on the pressure distribution for any of the three angles of attack at M = 6.8; this 
suggests attached flow on the lower surface. The blast -wave effect combines with the 
power law curvature of the lower surface to cause relatively constant pressure over the 
forward two -thirds of the surface at -8° angle of attack, over the forward one -half of the 
surface at 0° angle of attack, and over the forward one -third of the surface at 15° angle 
of attack. To the rear of these constant-pressure areas, the pressure trend is toward 
increasingly severe adverse pressure gradients at higher angles of attack. 

The heat-transfer results at M = 6.8 are little affected by the varying Reynolds 
number on the part of the concave surface where the pressure is constant, but they are 
greatly affected by the changing Reynolds number in the region of adverse pressure gra- 
dient. For -8° and 0° angles of attack, the heat-transfer data at the lowest Reynolds 
number suggest laminar flow, and the data at increasing Reynolds numbers suggest tran- 
sitional flow. At 15° angle of attack, none of the data suggested laminar flow over the 
entire surface. Increasing Reynolds numbers were accompanied by an apparent forward 
movement of the start of transition. 

The theoretical heat transfer to the concave surface in the region of constant pres- 
sure is another example of variance between results obtained from the two theoretical 
methods. The disagreement is most apparent for -8° angle of attack. It is also obvious 
that the theory of reference 6 involving the machine calculation is in better agreement 
with both the level and the trend of the heat -transfer data. Remember that the integral 
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associated with the transformed streamwise coordinate in the theory of reference 6 is 
evaluated from the stagnation point of the blunt leading edge by the use of the Newtonian 
pressure distribution around the cylinder, whereas the closed-form solution by the pres- 
ent author (ref. 3) does not account for this leading-edge bluntness in the analogous inte- 
gral. In this instance the failure of the assumed power law pressure distribution to cor- 
rectly predict pressures in the leading-edge region causes inaccuracy in this theoretical 
heat-transfer calculation. At the higher angles of attack of 0° and 15°, the two theoret- 
ical heat -transfer calculations are observed to approach each other at increasing angles 
of attack. This is expected since the integral of a higher nearly constant surface pres- 
sure function is less affected by the relatively smaller pressure perturbation near the 
leading edge. 


CONCLUDING REMARKS 


Pressure and heat transfer have been measured on convex and concave curved sur- 
faces following a hemicylindrically blunted leading edge at Mach numbers of 6.8 and 9.6 
and Reynolds numbers based on length from 0.5 x 10® to 3.0 x 10®. The curvature of 
both surfaces was specified by the same power law relationship. 

Power law functions were determined which fit the pressure data. The empirical 
constants were then used to calculate the heat transfer by a previously published closed- 
form solution by the present author (AIAA Journal, April 1967). This method used 
Prandtl number 1 (Np r = l) solutions of the laminar boundary layer with an arbitrary 
pressure gradient. These solutions were reduced to a closed form by assuming that the 
pressure varies as the power of the sum of the streamwise coordinate and an arbitrary 
constant or p (x + A) n . This condition permitted the evaluation of the streamwise 
coordinate transformation. Then the Np r = 1 solutions are expressible in a simple 
closed form in terms of physical quantities in the real plane. These theoretical values 
were compared with a point -by -point solution of the laminar boundary -layer equations 
by an iterative electronic computer method with similarity assumed and with the experi- 
mentally determined heat transfer. In most cases the two theories and the data were in 
good agreement. In some cases the simple theory of the present author gave less than 
desirable accuracy and certainty. The failure of the simple theory in these instances 
was attributed to the failure of the assumed power law to predict the correct pressure 
in the region of the leading edge. An example was shown where the assumed power law 
which adequately agreed with pressure aft of the leading edge was perturbed slightly with 
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a constraint to predict the stagnation line pressure accurately, and where the accuracy 
of the calculated results was improved. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 13, 1967, 

129-01-07-07-23. 


16 



APPENDIX A 


A SUMMARY OF THE THEORETICAL RESULTS OF REFERENCE 3 

The details of a simple theoretical method for calculating characteristics of the 
laminar boundary layer are shown in reference 3. This method involves the assumption 
of p cc (x + A) n , boundary -layer growth in a constant entropy field, and a perfect gas 
equation of state. These assumptions are applied to the Np r = 1 compressible solu- 
tions of the laminar boundary layer, that are shown in reference 4. 

The boundary -layer thickness and the skin friction are little affected by the 
Npj. = 1 assumption since the solution of the energy equations has little effect on these 
parameters. The heat-transfer calculations are directly involved with the solution of 
the energy equation and are corrected for Np r = 1 by the method suggested in refer- 
ence 5. The final relations which involve no hypersonic assumptions appear as follows: 

The relationship between the Falkner-Skan and the Li-Nagamatsu pressure gradient 
parameters is 


m = 


n(l - r) / I - g \ _ n(l - y) , 

2y 
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The local and total skin-friction coefficients presented in a more accurate form 
than that seen in reference 3 are 



The values of Kj to K 5 may be evaluated in the frame of reference 2 without hyper- 
sonic approximations (with Kj, K 2 , and K 4 appearing in more accurate form) and 
appear as follows: 
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The relations for 9 and K5 were not defined or discussed in reference 2 or 3 
and are therefore discussed as follows. 

In a manner analogous to that given for the displacement thickness in reference 3, 
the expression for the momentum thickness is as follows: 


9 = 



(A16) 


If constant pressure is assumed across the boundary layer, the expression for the 
momentum thickness in a perfect gas for Np r = 1 becomes 


9 = 



f’(l - f')d?7 


(A17) 


Use of definitions from reference 2 allows equation (A17) to be simplified as follows: 





v tr 


(A18) 


or 
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(A19) 


If the assumption p cc (x + A) n is made and the procedure shown in reference 3 is fol- 
lowed, the momentum thickness may be expressed in terms of real quantities 



(A20) 


When this equation is cast into the frame of reference 2, 




(A21) 
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where 9 represents the zero pressure gradient value with local conditions the same as 
the free stream. Then 
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APPENDIX B 


THEORETICAL PRESSURES ON BLUNT CURVED PLATE 


Numerous methods have been advanced for the determination of pressures on sharp 
or blunt plates in hypersonic flow. A method suggested in reference 14 has been tried 
for the calculation of the surface pressure in the present case. The method used was the 
linear combination of the local wedge pressure with the contributions expected for the 
blast -wave effect from the blunt leading edge and for the viscous effect from the growth 
of the displacement thickness as follows: 


p _ Pw + ^Pb + ^Py 

Poo " Poo + Poo + Poo 


(Bl) 


The local wedge pressure was determined by the straightforward use of two- 
dimensional oblique shock or Prandtl-Meyer expansion theory for the pressure on a 
sharp flat plate at each local surface angle. The bluntness contribution was calculated 
by using equation (7) of reference 13 together with a correction suggested in reference 15 


4p b 

p» 
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0.187e|^7(y - 1)- 


M*C 


x/d 


D,Ze | 


2/3 


0.26 


(B2) 


The viscous contribution was evaluated by the use of what was referred to as the "com- 
plete theory" in reference 15. This involved the use of equation (15) of reference 15 
which is (in the terminology of this report) 



r(r - 1 ) 

4 


K 2 5 + yK 5 




(B3) 


With equation (7) of reference 15, which is (in the terminology of this report) 

K§ = -sp -p====-(l - n ) (B4) 

where K4 was evaluated as suggested in reference 15. The calculations showed the 
viscous contributions to be negligible except at low pressure ratios. 

The comparisons of these theoretical values with measured pressures are shown in 
figure 11. Comparisons on the convex surface are seen in figures 11(a) and 11(b) and 
agreement of level and trend is good. It is apparent, however, that agreement between 
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theory and data near the leading edge is better for a = 15° than for the other two angles 
of attack. For a = 15°, the surface is at about zero inclination near the leading edge. 

Comparisons of data and theory on the concave surface at M = 6.8 are shown in 
figure 11(c). In this figure, the linear combination of the two-dimensional expansion 
theory with the blast wave and viscous contribution has been less successful than for the 
convex surface. The poorest agreement is in the neighborhood of the nose for the wing 
at a = 0° and a = -8°. For these angles of attack, a blunt leading edge is followed by 
a surface at a negative angle of attack not near zero. This failure of this method is 
pointed out in reference 14. 

The pressure measurements on the concave surface at M = 9.6 showed the same 
trend as those at M = 6.8 but had greater scatter and are not included in figure 11. 
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L Temperature and pressure leads 



(a) Detail of concave surface. 

Figure 4.- Heat transfer and pressure model. 
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(b) Rj = 7.3 x 10 6 . 
Figure 6.- Concluded. 
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(a) Convex surface. 


Figure 7.- Oil-flow photographs at M = 6.8 and R[ - 7.3 * 10 6 . 
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(c) Vortex evidence alongside of model. L-67-6611 

Figure 7.- Concluded. 
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(a) Convex surface; a = -8°. 

Figure 8.- Pressure and heat transfer on blunt curved plates with and without end plates at M = 6.8 and Rj 
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(b) Convex surface-, a = 0°. 


Figure 8.- Continued. 




(c) Convex surface; a = 15' 
Figure 8.- Continued. 
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(e) Concave surface; a = 0°. 
Figure 8.- Continued- 
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(f) Concave surface; a = 15°. 
Figure 8.- Concluded. 
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(d) a = 0°; M = 9.6. 
Figure 9.- Continued. 
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(b) a = 0°; M = 6.8. 
Figure 10.- Continued. 
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(d) a = 15°; M = 9.6. 
Figure 10.- Concluded. 





(a) Convex surface; M = 6.8. 

Figure 11.- Comparison of simple hypersonic approximations with measured pressure distributions 
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(c) Concave surface; M = 6.8. 
Figure 1L- Concluded. 
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